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Abstract. We prove that the density fluctuations for a zero-range process 
evolving on the supercritical percolation cluster are given by a generalized 
Ornstein-Uhlenbeck process in the space of distributions S'(R d ). 



1. Introduction 

Consider the infinite cluster of a supercritical bond percolation model. On this 
random graph, we define a zero-range process, which can be defined as a system 
of symmetric, simple random walks on which the hopping time of particles at a 
given site depends only on the number of particles at that site. It has been recently 
proved that the scaling limit of a simple random walk on the percolation cluster is 
given by a Brownian motion with a diffusion coefficient that does not depend on the 
particular realization of the percolation cluster [lOj , [3] . With this result in mind, it 
is natural to raise the question about the collective behavior of a system of random 
walks evolving on the percolation cluster, possibly with some interaction. More 
precisely, we want to study the density fluctuations of this model. Scale the lattice 
by 1/n and give a mass 1/n to each particle. In this way we obtain a measure in 
R d , which we call the empirical density of particles. Under a diffusive time scaling, 
it has been proved for a simple exclusion process [5] that the empirical measure 
converges in a proper sense, to a solution of the heat equation dtu — DAu, where 
D is the diffusion coefficient of the underlying random walk. 

This result can be interpreted as a law of large numbers for the density of par- 
ticles, and therefore the question about the central limit theorem arises. For the 
simple exclusion process, duality techniques allow to reduce the law of large num- 
bers and the central limit theorem for the density, to suitable problems for simple 
random walks. In order to treat a more general case, in these notes we discuss 
the density fluctuations for the zero-range process. Due to the inhomogeneity in- 
troduced by the randomness of the percolation cluster, the model turns out to 
be non-gradient. Following our previous work [6], a functional transformation of 
the density empirical allows to put the model back into the setup of gradient sys- 
tems. The main tool allowing to make this functional transformation is a form of 
the compensated compactness lemma of Tartar, which is obtained in [5 using two- 
scale convergence. Therefore, the proof of a central limit theorem for the density 
fluctuation is reduced to the so-called Boltzmann-Gibbs principle, which roughly 
states that fluctuations of non-conserved quantities are faster than fluctuations of 
conserved quantities. In the right scaling limit, only fluctuations of the density of 
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particles survives, allowing to obtain a martingale characterization of the scaling 
limit. In this way we prove that the limiting density field is given by a generalized 
Ornstein-Uhlenbeck process, as predicted by the fluctuation-dissipation principle. 
Unfortunately, and only in dimension d = 2, a regularity result for the solutions of 
the heat equation on the percolation cluster prevent us to obtain such a result. 

The main technical innovation in this paper is what we call the connectivity 
lemma. In [6] , an ellipticity condition allows to compare the relaxation properties 
of our system with the relaxation properties of a zero-range process in the absence of 
the random environment. But the percolation cluster do not satisfy this ellipticity 
condition, since by construction the jump rate through a closed bond is 0. Moreover, 
the percolation cluster is locally non connected, since for a typical box of fixed 
size, the part of the percolation cluster lying inside this box is not connected. 
The connectivity lemma says that enlarging a bit the box we obtain a connected 
graph with good ergodic properties, allowing to get bounds on the relaxation to 
equilibrium. 

We point out that our results are still true in a more general context. We have 
chosen the zero-range process on the percolation cluster for simplicity and to capture 
the essential difficulties. Take for example a random barrier model on X d on which 
the conductances are not bounded below (see [5] for a more detailed definition), 
and take an interacting particle system satisfying the gradient condition (before 
randomizing the lattice). Our results applies as well for these models, leading to 
the same results. 

These notes are organized as follows. In Section [2] we introduce the model and 
we state the main results. In Section [3] we introduce the corrected density field 
and we show how to obtain the main results from the compensated compactness 
lemma and the Boltzmann-Gibbs principle. In Section [4] we state and prove the 
connectivity lemma and we obtain the Boltzmann-Gibbs principle starting from it. 

2. Definitions and main results 

2.1. The supercritical percolation cluster. Let £ = {e = {xy);x,y £ Z d , \x — 
y\ = 1} be the set of nearest-neighbor, non oriented bonds in Z d . Let u> — {w(e); e <E 
£} be a sequence of i.i.d. random variables with P(w(e) — 1) = 1— P(u>(e) — 0) = p. 
Let us fix a realization of the sequence u>. Whenever uo(e) = 1, we say that the bond 
e is open. Otherwise we say that the bond e is closed. We say that two sites x, y in 
Wi d are connected, which we denote by x <-> y, if there is a finite sequence of points 
{xq = x, #1, . . . , x n = y) such that for any i £ {0, . . . , n — 1}, |a;<+i — Xi\ = 1 and 
the bond (xiXi+i) is open. For each x £ T, d , we define C(ui,x) ={t;e Z d ;x <-> y}. 
We call C(uj,x) the cluster containing x. It is well known [7] that the probability of 
C(u>, x) being infinite is an increasing function of p, and that there exists p c £ (0, 1) 
such that this probability is strictly positive for p > p c and equal to for p < p c . 
It is also well known that for p > p c , with probability one, there exists a unique 
cluster of infinite cardinality. We denote this cluster by C(u). We call C(u>) the 
supercritical percolation cluster. From now on, we fix a number p £ (p c , 1] and a 
sequence u> for which C(u>) is well defined. We also define 9(p) = P(0 £ C(ui)). 

2.2. The zero-range process in C(cu). For x, y in Z d , we say that x ~ y if 
\x — y\ = 1 and (xy) is open. Let g : No = {0, 1, ...} — * [0, oo) be a function with 
g(0) = 0. The zero-range process in C(u>) with interaction rate g(-) is the Markov 
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process £ t in fl = N generated by the operator L given by 

Lf{o= ff(^))[/(r t, )-/(o], 

x,yel. d 

where / : il — > R is a Zoca/ function, that is, it depends on £(ai) only for a finite 
number of sites x G C(u>), and £ x ' y is given by 

!£(x) — 1, z = cc 
£(z), z^a;,y. 

The dynamics of this process is easy to understand. At each time, a particle 
jumps from a site x G C(u>) to a neighbor y ~ x with exponential rate g(k), where 
k is the number of particles at site x at that time. When the initial number of 
particles is finite, the process £t is a continuous-time Markov chain in O. For the 
construction of the process in the case of an infinite number of particles, we refer 
to ffl. 

In order to have a well defined family of ergodic, invariant measures for the 
process £t, we assume that gin) > for any n > 0, and that 

sup \g(n + 1) - g(n)\ < +oo. (2.1) 

n 

Under these conditions, g(-) is bounded by a linear function: there exists cq such 
that g(n) < c$n for any n. Let us define the uniform, product probability measures 
in f2 by the relation 

1 ^ 
Z(ip) g(k)\ 

where g(k)\ = g(l) ■ ■ ■ g(k) for k > 1, g>(0)! = 0, i/3 G [0, oo) and Z(<p) is the 
normalization constant. Due to l|2.ip . these measures are well defined up to some 
critical value (f c . Notice that the parameter ip has an interpretation as the strenght 
of the interaction, since J g{^{x))i' v (d^) — ip. It is not hard to prove that the 
measures {v v ;tp G [0, y c )} are ergodic and invariant under the evolution of £t. 
Since we are interested in density fluctuations, we will reparametrize this family of 
measures by its density of particles. For that purpose, define p[}p) = J ^(x)D ip (d^). 
It is not hard to see that tp > p{f) is a diffeomorphism from [0,ip c ) to [0, oo). In 
particular, the inverse function p i— * <p{p) is well defined. We define v p = Vm( p ), 
measure for which now J S,(x)v p (d^) = p. 

2.3. The density fluctuations. Let p G (0, oo) be fixed and consider the process 
£™ = £tn 2 starting from the measure v p . We denote by P n the distribution of the 
process £™ in the Skorohod space 2?([0, oo), f2) of cadlag trajectories in fi, and by 
E n the expectation with respect to P n . 

Denote by C c (R d ) the set of continuous functions G : K. d — > M with compact 
support. By the ergodic theorem, it is not hard to see that 

lim — d V G{x/n)i{x) = 6(p)p f G{x)dx, 

almost surely with respect to the probability measure P ® v p . In the previous 
result, the ergodic theorem is invoked twice: first to state that the density of points 
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belonging to C(u>) when properly rescaled is equal to 9{p), and then to state that 
the density of particles is equal to p. Since the measure v p is invariant under the 
dynamics of £, the same result is also valid if we replace £ by £™ in the previous 
expression. We are interested on a version of the central limit theorem for this 
quantity. Let us define, for each test function G, the density fluctuation field y™ (G) 

by 

For topological reasons, it will be convenient to restrict the previous definition 
to functions G G iS(M d ), the Schwartz space of test functions, although y™(G) 
makes sense for more general test functions. The process y™ defined in this way 
corresponds to a process on the Skorohod space T>([0, oo), <S'(M d )), where S'(R d ) is 
the space of distributions on M. d . Now we are ready to state our main result: 

Theorem 2.1. Fix a particle density p > 0. In dimension d > 2, for almost 
all lu, the sequence of processes {^"In converges in the sense of finite- dimensional 
distributions to a generalized Ornstein- Uhlenbeck process yt of mean zero and char- 
acteristics ip'(p)DA, y' D(p(p)V , where D is the limiting variance of a symmetric 
random walk in C(u>). 

Notice that this theorem implies in particular a central limit theorem for the 
sequence {y?{G)} for any G G S(R d ). 

3. The corrected fluctuation field 

3.1. The corrected fluctuation field. The proof of Theorem 12.11 is based on 
Holley-Stroock's characterization of generalized Ornstein-Uhlenbeck processes [5]. 
For each t > 0, let T t be the cx-algebra on 2?([0, oo), <S'(IR d )) generated by the 
projections {y s (H);s <t,H e S(R d )}. The process yt admits the following char- 
acterization: 

Proposition 3.1. There exists a unique process y t inC([0,oo),5'(R' 1 )) such that: 

i) For every function G G S(R d ), 

M t {G) = y t (G) - y (G) - <p'(p)D f y s (AG)ds 

Jo 

and 

{Mt{G)) 2 — 6(p)(p(p)Dt [ (VG{x)fdx 

are T \ -martingales. 

ii) 3^o *s Q Gaussian field of mean zero and covariance given by 

E[y (G)y (H)] = 9(p) x (p) f G{x)H(x)dx, 

where X (p) = Var(?/(0); t^) and G,H G S(R d ). The process yt is called 
the generalized Ornstein- Uhlenbeck process of mean-zero and characteristics 
ip'(p)DA and ^ 'ip(p)DV '. 

Now the idea behind the proof of Theorem 12.11 is simple. We will prove that 
the sequence of processes {y™} n is tight and that every limit point satisfies the 
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martingale problem stated in Proposition 13.11 With these two elements in hand, 
we will be able to conclude Theorcm l2.ll By Dynkin's formula, for each G E 6>(R d ), 

M?(G) = X l (G) - yS(G) - [ n 2 Ly?(G)ds 

Jo 

is a martingale with respect to = er{£";s < t}. The quadratic variation of 
M t n (G) is given by 

C n 2 {Ly™(G) 2 - 2y;(G)Ly:(G)}ds. 

Jo 

Notice that the second step (that is, that the limit points satisfy the martingale 
problem) requires to replace, in some sense, n 2 Ly™(G) by ip'(p)Dy™(AG), and 
to replace Ly™{G) 2 - 2yi l (G)Ly™(G) by 9(p)ip{p)D J Rd (VG{x)) 2 dx. Let us take 
a more careful look at these two terms. We start with the second one. Simple 
computations show that 

LX l (G) 2 - 2y-(G)Ly-(G) = i E g(€(x))n 2 (G(y/n) - G{x/n))\ 

xec{uj) 

y:x~y 

By the ergodic theorem, when integrated in time this quantity should converge to 
9(p)ip(p)t J (VG(i)) dx. Notice that we have missed the factor D in the previous 
computation. We will return to that point later. Now let us take a look at the first 
term: 

^"( G ) = ^ E (s(CW)-vW)4G(i/«), 

£CGC(w) 

where g(-) is the interaction rate and C n is the generator of the associated symmetric 
random walk in C(lu). For a function F : R d — > R and x E C(uj), C n F(x/n) is defined 

by 

£ n F(x/n) = n 2 ^ F(y/ri) - F(x/ri). 

yeC(uj) 
yr^x 

The so-called Boltzmann-Gibbs principle will allow us to replace, when integrated 
in time, the expression g(££(x)) — ip(p) by ip'(p)^(x) — p) in the sum above, 
allowing us to write n 2 Ly^{G) as ip' (p)y t (£ n G) plus a term that is negligible as 
n — » oo. Now we can see what the problem is. Assume that {y™} n is tight. 
Take a limit point y^° of {3^}w We want to say that y^{C n G) converges to 
y^°(DAG) along the corresponding subsequence. But yf converges to 3^°° only 
on a weak sense. Therefore, we should need strong convergence of C n G, which is 
easily checked not to hold. 

The way to overcome this problem is to use the compensated compactness lemma. 
The idea is to choose, for each n, a test function G n in such a way that C n G n 
converges strongly to DAG as n — > oo. We therefore define the corrected fluctuation 
field as in [6]. Fix some A > 0. For each G E 5(R d ), define G n : C(w) -> R d as the 
solution of the resolvent equation 

XG n {x) - £nG n (x) = XG{x/n) - DAG(x/n). (3.1) 

Then the corrected fluctuation field y™' is given by 

xec(ui) 
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Notice that C n G n = A(G„ — G) + DAG. In particular, strong convergence 
of C n G n to DAG follows from strong convergence of G n to G. The following 
proposition tells us that this is, indeed, the case. 

Proposition 3.2 (Faggionato [S]). There is a set of P -total probability such that 
for any G G the sequence {G n }„ converges to G n in the following strong 

sense: 

lim 4? y \GJx) -G(x/n)\ 2 = 0. 

n—>oo n ^— ' 1 1 
x£C(uj) 

In particular, since the invariant measure v p is of product form, we conclude that 
y^ X (G n ) - y?(G) vanishes in £ 2 (P n ) as n — > oo. We will see that the standard 
scheme, tightness plus uniquenes of limit points via martingale characterization, 
can be accomplished for y™' X . 



3.2. The martingale problem for the corrected fluctuation field. When 
considering the c< 
resentation gives 

Ml 

In this formula, we have defined 



considering the corrected fluctuation field y"' instead of y?, the martingale rep- 
i gives 

C< A (G) = y?'\G) - y^\G) - [ 6?(A(G„ - G) + DAG)ds. 

Jo 



n d/2 

The quadratic variation of A/"' A is given by 



(M t "' A ) 



f ~ d Y. n 2 g(C(x))(G n (y)-G n (x)) 2 ds. 



y:x~y 



Multiplying the resolvent equation (13. ip by G n (x) and summing over x G C(u>) 
we see that 

n Um^ E n 2 (G n (y)~G n (x)) 2 = 6{p)D J (WG(x)) 2 dx. (3.2) 



xec(u) 

y-x~y 



Let us denote by Q x the distribution of the process y?' in V([Q, oo), S'(R d )). 
Following [6], it is not hard to prove that the sequence {Q„} is tight. In the same 
way, for d > 3, we can prove that 



lim (M t n ' A ) = <p(p)0(p)D / NG{x)) dx. 

n^oc J 

The arguments in [B], and only in dimension d = 2, require some regularity of 
G„ that is missing in our situation. This is the only point on this article where we 
require d > 3, otherwise the proofs do not depend on the dimension. Now we will 
state in a more precise way the Boltzmann-Gibbs principle. 

Theorem 3.3 (Boltzmann-Gibbs principle). For any t > and any F G C c (R d ), 



lim E r , 

n — >oo 



{e-(F)-^(p)yi\F)}ds 



DENSITY FLUCTUATIONS ON PERCOLATION CLUSTERS 



7 



As in [6], from the tightness of {Qn}, Theorem 13.31 an the expression (13. 2| allow 
us conclude the following result: 

Theorem 3.4. Fix a particle density p > 0. In dimension d > 2, for almost any 
ui, the sequence of processes {y™' X } n converges in distribution with respect to the 
Skorohod topology in X>([0, oo),<S'(R )) to a generalized Ornstein- Uhlenbeck process 
y t of mean zero and characteristics (p'(p)DA, Dtp(p)'S7 , where D is the limiting 
variance of a symmetric random walk in C(w). 

Theorem 12. II is an immediate consequence of this Theorem and Proposition [X2] 



4. The Boltzmann-Gibbs principle 

In this section we prove Theorem l3.3l The main point that makes the proof of this 
theorem different from the proof in [6] is the lack of ellipticity: by construction, the 
jump rate between two neighboring sites x, y is equal to when the bond (xy) is not 
open. In [5] , due to the ellipticity condition we could compare the generator of the 
process with the generator of a zero-range process without the random environment, 
but with a slower jump rate. Since the Boltzmann-Gibbs principle holds for the 
latter process, it should hold for the former, since the dynamics is faster. In our 
case, we do not have a proper slower process to compare with. 

Let us explain better the intuition behind Theorem 13.31 The idea is that non- 
conserved quantities fluctuate faster than conserved ones. Since the only conserved 
quantity for the zero-range process is the number of particles, it is reasonable that 
at the right scale, the only part of the fluctuation field 8™ that is seen at a macro- 
scopic level is its projection over the conservative field y™. We can think that 
non-conserved quantities equilibrate locally, while conserved quantities need to be 
transported in order to equilibrate. In particular, we will see that the form of the 
graph is not really important in Theorem 13.31 What is really important is the 
connectivity of the graph: if a graph has more than one connected component, 
then there is more that one conserved quantity: the number of particles on each 
connected component. 

For n S No, define V(n) — g(n) — Lp(p) — ip'(p)(n — p). The strategy of the proof 
of Theorem 13.31 introduced by Chang [4] is the following. Fix a positive integer k. 
Divide the support of the test function F into small boxes of size k/n. Since F is 
continuous, we can average the function V(S,(x)) over the corresponding boxes on 
the lattice. Then we use some sort of ergodic theorem to reduce the sum over many 
blocks integrated in time, into a sum over a single block. This last problem is a 
static one. A new ergodicity argument, now with respect to the invariant measure 
v p will allow us to conclude. 

An important property of the lattice is that it is locally connected, in the sense 
that the restriction of the lattice to any box is still a connected graph. This is no 
longer true for the percolation cluster: the restriction of the percolation cluster to 
a box is, in general, not connected. Although most of it belongs to a big single 
component, the rest is spread over many small connected components. Of course, 
all these connected components are connected by paths that pass outside the initial 
box. The point is that these paths can be chosen is such a way that they do not go 
too far from the original box. We will develop these ideas in the next paragraphs. 
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4.1. A connectivity lemma. In this section we will state a result that we call 
the connectivity lemma and we will prove Theorem 13. 31 starting from it. The proof 
of this connectivity lemma is postponed to the next section. 

Let k, I be fixed positive integers. We will send k and I to infinity after n. 
We introduce two intermediate scales in our problem as follows. For simplicity we 
assume n/k 6 N, and pasting a sufficient number of cubes, we can assume that 
the support of F is contained in the cube (S, 1 — S) d for some <5 > 0. We can 
split the cube A„ = {1,2,..., n} d into (n/k) d non-overlapping cubes of side k. Let 
{Bj,j — 1, . . . , (n/k) d } be an enumeration of those cubes. Define Bj as the box of 
side (21 + l)fc, centered at Bj. In particular, the cubes Bj are the union of (21 + l) d 
cubes in {B®}i, except for the ones near the border of A n . Now define 

Bj = BjnC(w), 

Bj = {x E Bj-.x <-> B°}. 

In other words, Bj is the portion of the cluster C(ui) inside the box Bj, and 
Bj is the portion of the cluster C(oj) inside Bj and connected to the box B ( j. We 
say that a cube Bj is good if Bj is connected. In other words, Bj is good if the 
connected components of B® are connected between them by paths that lie entirely 
in Bj. We will denote by B n the union of good cubes. The numbers k and I will 
be fixed most of the time. Therefore, in order to keep notation simple we do not 
make explicit the dependence of B n in k and I. A cube that is not in B n will be 
called bad. We will call Bj indistinctely the set of points already defined, and the 
subgraph of C(oj) corresponding to these points. 

Define C n (ui) = C(uj)PiA n . Theorem [33] is a consequence of the following lemma: 

Lemma 4.1 (Connectivity lemma). For each e > there exists I > such that 

i) lim limsupE„|Y f —L V V(£(x))F(x/n)ds) 2 | =0 (4.1) 

••M- r \B n \ 

it) hm sup lim sup -5 < e, (4.2) 

where \C n (uj) \ B n \ denotes the cardinality of the set C„(lo) \ B n . 

The first part of the lemma says, roughly speaking, that Theorem 13.31 is true if 
the considered graph is locally connected. The second part says that C(uj) is not 
locally connected only on a small portion of the lattice. 

Now let us prove Theorem 13 . 31 assuming Lemma |4. II The expectation appearing 
in the statement of Theorem 13.31 is bounded by 



2E„ 



' ^ V(C(x))F(x/n)ds' 2 



n d/2 

xeB n 



2E„ 



? d/2 



£ V(C(x))F(x/n)ds 



2i 



The first expectation goes to as n — > 00 and then k — > 00 by Lemma |4.1[ part 
i). By Schwartz inequality, the second expectation is bounded by 



t 



2 

x£C n (uj)\B n 



F(x/nfE p [V(i(x)f 
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which turns out to be bounded by C(t, F, p)\C n {ui) \ B n \/n d . This last expression 
vanishes as n — > oo and then k — > oo by Lemma 14. 11 part n), which proves Theorem 



4.2. Proof of the connectivity lemma, part i). In order to simplify the nota- 
tion, we also denote by B n the set of indices j in {1, . . . , (n/k) d } for which B® is a 
good cube. Let us start with part i). We will manipulate the term 

f J/2 E VKW)f(^)*) 2 ] (4-3) 

until we arrive to an expression that does not depend on t. Taking the positive 
and negative parts of F, we can assume, without loss of generality, that F is non- 
negative. For each j, take a point yj in Bj. Since the function F is uniformly 
continuous, we can rewrite the integrand in (|4.3p as 



^ E F ^ln) £ V(ff(*)) (4.4) 

plus a rest that vanishes in C 2 (v p ) as n — > oo and then fc — » oo. Notice that now 
we have introduced an averaging of the function V(£(a;)) over boxes of side k. 

Now we explain the point where we make use of the time average in (|4.3[) . For 
each j, denote by L B . the restriction of the generator L to the set Bj. Observe that 

the zero-range process restricted to Nq j is ergodic on the set of configurations with 
a fixed number of particles, exactly due to the fact that Bj is a good cube. For any 
two functions /, h : fi — ► R, denote by (/, h) p the inner product with respect to v p . 
We define by 

||/||?= (f,-Lf) p , 
\\f\\ 2 _ 1 ^snp{2(f,h) p -\\h\\l} 

h 

the Sobolev norms associated to L. Here the supremum is taken over functions 
h G £?{v p ). Take an arbitrary (by now) family of functions {fj',j G B n } with 
fj : Nq 3 — > M. We have the following Sobolev inequality (Prop. Al.6.1 of [S]): 



E„ 



fiE ^n)Wi(e)^) 2 ] < ^11-4 E ^(%/^)^/iH 2 -r 

Let us call W n the term inside the norm. We will bound HW^H^ using the 
variational formula of 1 1 • | \ 2 _ 1 introduced above. Then, we need to estimate (h, W„) p 
for h given. By the weighted Schwartz inequality, 

J hL Bj fdv p < l-(f,-L Bj f) p + ^-(h,-L Bj h) p . 

Choose 7i = n 2+d / 2 /F(yj/n) and observe that J2j(h, —Lsjh) p < (h,—Lh) p . 
Plug this estimate into the variational formula for 1 1 Wn 1 1 ?_ x to discover that 

E «[( fiE F(y 3 /n)L B]fj (C)ds) 2 ] < ^ F(y 3 /nf(f 3 ,-L B] f) p . 

Notice that, for fj fixed, the right-hand side of this inequality is of order 1/n 2 . 
Notice also that we can not take the same function fj for different j, since the 
clusters Bj are different for each j. However, k and / are fixed. In particular, the 
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possibilities for the cluster Bj are many, but finite. Let S be the set of connected 
graphs contained in a box of size (21 + l)k. Let us denote by 77 a generic element 
(that is, a graph) in S. Take now a family of functions {/ 7? ; 77 G S} and define 
fj = p if Bj = rj. With this notation, we can bound the right-hand side of the 
previous inequality by 



40*11*11; 

„d+2 



iie-E 



where q n (r)) is the number of boxes in B n with graph 77, and is the generator L 
restricted to the graph 77. Since q n (r])/n d is bounded, the functions {f n } n are fixed 
and also S is fixed, this last quantity goes to as n — > 00. In particular, this means 
that we can discount a sum of the form n~ d / 2 J^. F(yj/n)LB j fj in (|4.4[) . We can 
also take the infimum over the families {f v } n , although only after n — > 00. In other 
words, we have reduced (|4.3p to the verification of 



lim sup lim sup E r , 

k~*oc [fny n— »oo 



1 



The good point is that in this sum the time integral has already played its part, 
allowing the introduction of the functions fj . Bound the expectation using Schwartz 
inequality twice, once to get rid of the time integral, and once more to make use 
of the product form of the measure v p . Notice that the boxes Bj at a distance 



greater than (21 - 
is bounded by 



l)fc are independent. Then, the expectation in the previous limit 



t 2 (2l + l) d \\F\ 



Remark that the sum above depends on n only through the number q n (rf)/n d . 
By the ergodic theorem, q n (rf) / (n / k) d converges, as n — > 00, to q(rj), which is the 
probability of a given subgraph Bj of C(ui) be (good and) equal to rj. Therefore, 
we are left to prove that 



lim sup 



t 2 (2l + l) d \\F\\ i c 



0. 



(4.5) 



Here we remark that the whole construction we have done has as purpose to 
get a connected graph in the sum above. For each graph rj fixed, and due to the 
ergodicity of the zero-range process on sets with fixed number of particles, 



^^[(Ev^-iv 



= E n 



where f denotes the number of particles of the configuration £ on the graph rj. 
Define, for pf > 0, tp(p') = E p [v(t(x))]. Notice that ip(p) = 4>'(p) = 0. This fact 
combined with the equivalence of ensembles allow us to show that 



sup E p 



x£rj 



< 
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where the supremum is over k > and rj £ S. Therefore, we have proved that the 
expression in (14. 5| is bounded by 

c( *"3 g+ "' !>>■ 

which goes to as k — > oo, since the sum in 77 is bounded by 1 (remember that q{rf) 
is a probability). 

4.3. Proof of the connectivity lemma, part ii). We will use a result of [2J, 
which roughly states that the percolation cluster has good connectivity properties. 

For any two points x,y £ C(uj), we define the distance D(x,y) as the length of 
the minimal path connecting x and y: 

D(x, y) = inf{n; 3x = xo, ■ ■ ■ ,x n = y with (x i - 1 x i ) open}. 

We have the following result: 

Proposition 4.2 (Antal-Pisztora 2 ]). There exists a constant 7 = 7(p, d) 6 [1, 00) 

limsup — logP(0 <-» z, D(0, z) > 7|z|) < 0. 

|z|->oo Pi 

The previous Proposition is telling us that, with high probability, two points that 
belong to C(lo), both inside a box of size k, are connected through a path contained 
on the box of size (7 + <5)fc, centered on the original box. This means that for I > 7, 
(|4.2p should hold. Now we will make this point more precise. 

For simplicity, assume k odd. Fix e > 0. Let T° be the box of size k, centered at 
the origin. Define = T° n C(u>), the intersection of T° with the infinite cluster. 
By the ergodic theorem, |rfe|/fc" — > 6(p) a.s. and also in probability. Fix fco such 
that Fk 7^ with probability bigger that 1 — e/2. By Proposition ^. 21 there is 5 > 
such that 

P(0 ^ z, D(0, z) > j\z\) < exp{-5\z\} (4.6) 
for \z\ is big enough. Taking a bigger k if necessary, we can assume that this 
inequality holds for any z such that \z\ > k Q . Fix x £ T° and consider for a 
moment the measure P(-\x £ C(ui)), well defined since P{x £ C(uj)) — 0(p) > 0. 
Fix I > 0. In analogy with our previous definition, we say that the box Tfc is good 
if each connected component of Tfc is connected to x through a path contained in 
^(2/+i)fe- This i s tne case i 1 aim on ly ^ eacn point of C(w) n9r° is connected to x by 
a path contained in r° 2;+1 ^ fe . When we do not fix a point x £ , we say that 
is good if any two connected components of F^ are connected by a path contained 
in ^%i+i)k' ^ > eacn one °f the points in C(ui) D (9r° are at distance from 
x at least \k — ko\ and at most \k + ko\- By (|4.6|) . the probability (with respect 
to P(-|a; £ C(uj))) of all these points being connected to a; by a path contained in 
r (2/+i)fc is bigger than 

1 - 2dk d ~ 1 exp{-S\k - k \} (4.7) 
as soon as I > r y+(l+j)k a /k. This is the case if, for example, I > 37. The quantity 
in (|4.7p goes to as k — > 00. Since this estimate is independent of the choice of site 
x e r° o , we conclude that 

limsupP(r/£ is not good ) < e/2. 
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Notice that the definition of goodness of a given box is translation invariant. In 
particular, with the Bj's notation, 



limsup P(Bj is not good ) < e/2. 

By the ergodic theorem, 



k — >oo 



,. \C n {u)\B n \ , . 

hm 2 = P\Pj is not good) 

in probability and a.s. We conclude that for I > 37, 



r r \CnH\B n \ 

lira sup lira sup S e > 

k— >oo n-^oo TL 



which ends the proof of (14.2 
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